Accelerated Monotonic Algorithms for Transmission Tomography
نویسندگان
چکیده
We present a framework for designing fast and monotonic algorithms for transmission tomography penalizedlikelihood image reconstruction. The new algorithms are based on paraboloidal surrogate functions for the loglikelihood. Due to the form of the log-likelihood function, it is possible to find low curvature surrogate functions that guarantee monotonicity. Unlike previous methods, the proposed surrogate functions lead to monotonic algorithms even for the nonconvex loglikelihood that arises due to background events such as scatter and random coincidences. The gradient and the curvature of the likelihood terms are evaluated only once per iteration. Since the problem is simplified, the CPU time per iteration is less than that of current algorithms which directly minimize the objective, yet the convergence rate is comparable. The simplicity, monotonicity and speed of the new algorithms are quite attractive. The convergence rates of the algorithms are demonstrated using real PET transmission scans.
منابع مشابه
Accelerated Monotonic Algorithms for Transmission Tomography - Image Processing, 1998. ICIP 98. Proceedings. 1998 International Conference on
We present a framework for designing fast and monotonic algorithms for transmission tomography penalizedlikelihood image reconstruction. The new algorithms are based on paraboloidal surrogate functions for the loglikelihood. Due to the form of the log-likelihood function, it is possible to find low curvature surrogate functions that guarantee monotonicity. Unlike previous methods, the proposed ...
متن کاملFast Monotonic Algorithms for Transmission Tomography
We present a framework for designing fast and monotonic algorithms for transmission tomography penalized-likelihood image reconstruction. The new algorithms are based on paraboloidal surrogate functions for the log-likelihood. Due to the form of the log-likelihood function, it is possible to find low curvature surrogate functions that guarantee monotonicity. Unlike previous methods, the propose...
متن کاملAccelerated Optimization Algorithms for Statistical 3d X-ray Computed Tomography Image Reconstruction
ACCELERATED OPTIMIZATION ALGORITHMS FOR STATISTICAL 3D X-RAY COMPUTED TOMOGRAPHY IMAGE RECONSTRUCTION by Donghwan Kim Chair: Jeffrey A. Fessler X-ray computed tomography (CT) has been widely celebrated for its ability to visualize the anatomical information of patients, but has been criticized for high radiation exposure. Statistical image reconstruction algorithms in X-ray CT can provide impro...
متن کاملNon-monotonic Poisson Likelihood Maximization
This report summarizes the theory and some main applications of a new non-monotonic algorithm for maximizing a Poisson Likelihood, which for Positron Emission Tomography (PET) is equivalent to minimizing the associated Kullback-Leibler Divergence, and for Transmission Tomography is similar to maximizing the dual of a maximum entropy problem. We call our method non-monotonic maximum likelihood (...
متن کاملOrdered subsets algorithms for transmission tomography.
The ordered subsets EM (OSEM) algorithm has enjoyed considerable interest for emission image reconstruction due to its acceleration of the original EM algorithm and ease of programming. The transmission EM reconstruction algorithm converges very slowly and is not used in practice. In this paper, we introduce a simultaneous update algorithm called separable paraboloidal surrogates (SPS) that con...
متن کامل